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^ ■ Abstract 

, We consider a 5-D gravity plus a bulk scalar field, and with a 3-brane. The 

^ I Darboux transformation is used to construct some exact solutions. To do this we 

• reduce the system of equations, which describes the 5-D gravity and bulk scalar 

. field to the Schrodinger equation. The jump conditions at the branes lead to the 

I jump potential in the Schrodinger equation. Using the Darboux transformation 

with these jump conditions, we offer a new exact solution of the brane equations, 

CN ■ which represents a generalization of the Rundall-Sundrum solution. For simplicity, 

^ ' the main attention is focused on the case when Hubble root on the visible brane 

' is zero. However, the argument is given that our method is valid in more realistic 

^—1 . models with cosmological expansion. 

■ 

o ■ 
(N : 
o . 

■ 1 Introduction. 

Recently, E.E. Flanagan, S.-H.Tye and I. Wasserman [1] considered the 5-D gravity plus 
' a bulk scalar field with parallel 3-branes, for which the action is 

: 

X': S = J d'xdy^l (^^R - ^ - V{<f>)^ - ^ J^^ d'x^\a,i<P), (1) 

The system of units is chosen such that k = 1, where k is 5-D gravitational coupling 
constant. The coordinates are {x^,y), for < ;U < 3, the 6th brane is located at y = yb, 
Qab is the 5-D metric and g^^^^ is the induced metric on the 6th brane. The brane tensions 
(Tb and the potential V are the functions of the bulk scalar field = 4>{y). The authors of 
[1] seek solutions to the field equations of the form 

ds^ = dy^ + A{y) i-dt^ + e^^'dikdx'dx^) . (2) 



To solve the equations of motion, following the praxis in the supergravity theory, they 
define a superpotential by = W'^cj)), where the dot denotes the derivative in y 

and the prime - in 0, and H is the effective Hubble constant on the brane. If if = then 
one gets 

= l^'2/8 - W'^/Q. (3) 

Using this equation, the authors of [1] have suggested a simple way for obtaining 
solutions (see also [2, 3]). On the other hand, equation (3) implies a supersymmetry 
connection between V and W , so one can write down the Darboux transformation for 



these equations. This is because for one-dimensional systems ^ the supersymmetry and 
Darboux transformation are the same in a sense [4]. The aim of this work is to use the 
Darboux transformation to obtain exact solutions of the field equations, which can be 
derived by minimizing the action (1). For simplicity, the main attention is focused on the 
models without the cosmological expansion, although our approach also works in more 
realistic models with H ^ (see Sec.5). 

One should be quite careful carrying the Darboux transformation into the brane world. 
The reason is connected with the stick-slip nature of fields: the jump conditions at the 
branes are <p{yh) ~ ^t^iVb) ~ ^bi'^b)- So one must make sure that the Darboux transfor- 
mation be compatible with the jump condition. However, this can be done. Having done 
so, we have constructed a new exact solution of the brane equations, which is some sense 
generalizes the Rundall-Sundrum, henceforth RS solution. 

This paper is organized as follows. In section 2 we consider the Darboux transforma- 
tion and it's compatibility with the jump conditions on the brane. In the next section we 
present the new solution, which can be called the dressed RS-brane. The n — 2 Darboux 
iteration and the role of shape- invariant potentials are discussed in section 4. In section 
5 we show that the Darboux transformation can work in the case with the cosmological 
expansion and address some problems, characteristic of such models. Some concluding 
remarks are given in section 6. 



2 Darboux transformation and jump condition 

If if = and u = A/A then from (1) and (2) we get the following system ([!]), 

202 2^ . 2 2^ 
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-iiT.<my-y^), u'^\--v{<t^). (4) 
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The third equation has the form 

and can be obtained from the system (4), so it is enough to consider one. If we have a 
single brane, which is located at y = 0, then a{(f){y))5{y) = aQ5{y), where ctq = (t(0(O)) 
(we assume that 0(0) ^ 0). Introducing ip = ipiy) = A^, we get the Schrodinger equation 

i^^{v{y)-p5{y)-X)ij, (5) 

where 

v{y)^-lv{<l>{y)) + X, P = ^. 

Here A should not be taken for the cosmological constant (which is included into the 
potential T^(0)), being the usual spectral parameter, which we need in order to iterate 
the Darboux transformation (see below). 

What about of the inverse problem of restoration A{y), (j){y) and V^(0), knowing '?/'(t/), 
v{y) and A? Let consider this problem for large values of y, where the 5-function term in 
(5) is omitted. First of all, we assume that v{y) is such that 

/+0O 
dy {I + \y\)\v{y) < oo. (6) 
-oo 



1^ = (^(y), so this field is the solution of an ODE. 



The condition (6) is generally unnecessary for the use of the Darboux transformation, but 
in this paper we content ourselves with such potentials only. 

To reconstruct A{y) from ^(|/), one needs to ensure that ip{y) > for all y. The ground 
state wave function is the only wave function in the physical spectrum, possessing this 
property. Furthermore, it is the only solution of (5) (not only in the physical spectrum) 
which can be used to find the physical scalar field (f) for the potentials (6). The reason is 
as follows: if ip is the ground state wave function, then there exist such real that 

■?/'(?/) — > e~^^^ as y — > +oo, 

so from the first equation (4) one gets 

12 3« 3 

(f> = -— = --— log^(2/) ^ — > as y^+oo. 

On the other hand, one can use i/j as the solution of (5) with A < A^„^. Such a solution 
does not describe the bound state, but it can be positive for any values of y (one can 
choose the constants of integration such that this function does not have zeros, so A{y) 
can be constructed starting out of this ip), but in this case one gets 0^ < as |/ — >■ +oo. 

This is why throughout this paper we shall study only the ground states of (5). 

On the other hand, the Darboux transformation (henceforth DT) allows one to con- 
struct potentials with any bound state ad hoc. Therefore it is natural to consider the DT 
as a method for constructing the exact solutions of (4). This is the crucial point of this 
paper. 

Let ijji = ■ip{y;Xi) and ip2 = V^(2/;A2) are two solutions of (5). We call ipi the prop 
function for the DT which has the form [5]: 

, , (1) -02^'! - '01'02 (1) o / \ \ 

V'2^'02 - -, , V ^ > -2—\0gxlji, A2^A2. (7) 

ipi dy^ 

The DT (7) is an isospectral symmetry of (5): the dressed function Tp2^ is the solution of 
the dressed equation (5), namely 

= (^w A2) 

with the new (dressed) potential v'^^^ and the same eigenvalue A2^'' = A2. 
The transformation law for the prop function is ([4]) 

V'l - = ^ (^1 + C2 l%!{y'w) , = {v^'^ - P^'^Kv) - Ai) (8) 

with arbitrary constants Ci^2 (and, of course the same eigenvalue A^^-* = Ai). 

What about the ground state? If is the wave function of the ground state for the 
initial potential, and A2 is the next level in the physical spectrum, then after the DT 
(7) the dressed function ^2^^ will be the wave function of the new ground state [6]. In 
this case, the spectrum of the Hamiltonian with v^^^ can be obtained by the deleting the 
level Ai from the spectrum of the initial Hamiltonian (the case (i)). On the other hand, 
suppose ipi = ipi{y; Ai) is not really the wave function of the bound state of the initial 
Hamiltonian, but simply does not have zeros and ipi — > e"*"^ as 7/ — > ±00. Then after 
(7) one gets a new Hamiltonian with the new ground state with A = Ai and the wave 



function ip^^^ from the (8), with C2 — (the case (ii)). The latter alternative is used 
throughout this paper. 

The jump conditions at the brane are 



(9) 



where ipi±{y) are the solutions of the Schrodinger equation (5) with the same potential 
but different eigenvalues: 

A±^{v(y)-(35(y)-Xi)A±- (10) 
It is easy to see that conditions (9) are the same as the jump condition from the [1]: 

The following theorem shows that the DT preserve the jump conditions (9) 

Theorem. Let il>i{y) = {i^i+iy > ^),ipi-{y < 0)} with i = 1,2 be the solutions of (10), 
such that the jump conditions (9) are valid. If = j^/'i^., ^/'i. | is the prop function 

then the Darboux-dressed function ■^g'^^ — |^/'2+ 5^^2-1 (^^e prop function for ■02j_ is 
respectively) save the same jump conditions (9) with new "tension" P^^^ — —(5: 



41(0) = ^r_^(o) = ^^^^(0), ^^;^(o) - ^no) = +mo). 



,(1), 



■,(1) 



J,(l)/ 



;ii) 



The proof is a straightforward calculation, which we leave to the reader. As a hint, there 
are two useful formulas: 



fi±^v- Ai 



42 = (Ai - A2) V'. 



/(I) 



where ti± = V'i±/V'i±- 

Iterating this theorem, one can conclude that the n-fold DT results in dressed func- 
tions, which retain the jump condition (9) with P'^"'^ = (— l)"/?, provided that the n prop 
solutions {ipi±} {i = l,..,n) of the initial Schrodinger equation (10) and one dressing 
solution ipn+i± = ip of the same equation, with the same potentials but a different A, 
satisfy the jump condition (9). If this is the case, then one can use the Crum formulas [7] 
to present the result as 



D^n\y) 



vt\y) -2— log i^t\yl 



(12) 
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particular, one can see that the jump conditions from [1] are valid: 



(^A^±''^ , where ijj^^^ = d'^ip/dy^. In 



u^^\0) = (-1) 



n+l 



2an 



where u'^^ = A^]!^ /A^J^^ and (Tq is the tension on the initial brane. 

We shall use the case (ii) , for which the ground state of the potential v^^ is defined 
by the formula 

v-i:' = = (13) 

3 Single dressed RS brane. 

To use the DT (7), (8) we start with a simple initial potential v = fi^ = const (i.e. 
with V = 0, see (6)). This is the case of the Randall- Sundrum (RS) model for which 
V = -3/iV8- The solution of (5) with Ai = is 

^.^^a.^e'^y + b.^e-'^y, ^ Mv > 0), ^,_^My<o), (14) 

with the real positive constants a^_j_ > 0, > 0. 
Prom (9) one gets 



From (15) one can see that Qj^+ftj^ — a^_-|-6^_ and hence use the following representation: 

= sin^ CKj. , = cos^ . (16) 

Using (16) and (15) one gets the compatibility condition as 

— sin^Q;_^ < < cos^q;,^, (17) 

so \P\ < 2ii. 

Therefore, to construct new exact solutions via the DT, one needs to use (14) with 
the additional conditions (15) and (17). 

Now we choose (14) as the prop function. After the DT (7) we get 

where 

(1) 8/^^Qi±^i± 

with = v'^^\y > 0) and v^l^ — v^^\y < 0). This potential has the jump on the brane 

.«(0)-„-(0).-4,,(|-.^). 

which is zero if we choose the tension as 

(3 = '^^^-^7-^^ (18) 

If this is the case, then v^^ can be obtained from v^^ by transposing a^^^ and . 



We choose = l/'0i±,^ therefore as y — > 0: 

^(i)(0) - ^a)(o) = -/3(i)v^«(o) = +f3^\'\Q), 

as it should be, according to (11). 

Using ^[]^, one can calculate that the dressed metric has the form 



so A^^\y) e ''^/^/y^a^, asy ^ +00, aMdA'^^\y) e^'^l'^ j asy ^ -00. Tweaking 
(18), one gets a^, — so 



/l«(y) 



=-M|y|/2 



y — >• ±00. 



"1+ 



In this case, one can set = = to transform A^^^ into the exact RS solution of the 
system (4). We get the initial RS potential via /3 — > — /3. 
Let us return to the case of general position. Then 



, 2 ' 



and a simple calculations yields the well-known sine-Gordon potential V^^{(^±): 

4 



V^«(0±)-^¥cos 



x/3 



(0±(y) - 0o±) 



with 



0(y)± = 0o± + v^arctan [^-j,^^ 



where 0o± arc arbitrary constants. 



At last one can calculate the quantity (ctq^'') • The result is 

(^aj^)y = (/,J'(0) - (/)L^)(0) = //VSsin - «_) cos + «_) 
Note that in the case of tweaking (18) 



TT 



/3 = 2//cos2q;+, a- = q;+ + - (2A^ -M) , 



where A?" is integer. 



4 Two applications 

There arc many possible applications of these results. In this section we briefly consider 
the n = 2 dressing of the RS brane and the role of shape-invariant potentials. 



^By the way, we choose the constants of integration C\ = 1 and C2 = to be in the ground state. 



1. Double dressed RS model 

Another solution of the same Schrodinger equation with A2 = A*^ — has the form 

^^,±=a,^e'^^ + 6,^e-'^^ = V'lC?/ > 0), V'.- = < 0), (19) 

Using the theorem from section 2, one assumes that 



«=-=«,.+^^i^. 6,_=6„-^^±^. (20) 



Using (7) for (19), one gets 



= ^ , 



n± = (1/ - (a,±a,±e('^+-)^ - \J^,^e-^^^^')y) + (1/ + //) (a,±6,±e(''-'^)^ - a^J^.^e^^-^^^) . 

(21) 

The function (21) will be the prop function for the second DT. In order to create a new 
ground state with the "energy" A2 < 0, one should choose > > 0, a^^ > but 
ftji < 0. In this case ip^^J does not have any zeros (nor has it any singularities). Taking 
into account (20), we conclude that positive a^^ must be such that: 



-2- v2^ + /5/if/3>0 «2+ V2i^-/^/if/3<0 

If this is the case, then we can use (12) and (13) to proceed with the second DT, yielding 



One can verify that we have arrived in a brane with the right jump condition and tension 
a^^ — +(Jo. Note that the potential v^^^ has two levels with Ai = and A2 = /x^ — u"^. 

2. Shape-invariant potentials 

Another way to obtain exact soluble potentials via DT is to use the shape invariants 
[8]. If an initial potential is a function of y and some free parameters ai : v = v{y; ai), and 
after the DT one gets v^^^ — v^^\y] a-^^) then v is called a shape-invariant (SP) potential. 

SP-potentials are common in quantum mechanics. An example is the harmonic os- 
cillator [4]. We stress that exact soluble potentials from [1] (for the models without 
cosmological expansion) are SP-potentials too. We confine ourselves with considering 
three examples from the above article, retaining the terminology. 
A. Even superpotential. In the case of DFGK model [9] we get (with A{0) = 1): 

ay 



In this case 



After the DT 



logA(y) = -f + ^(l-e^^-). 



Q A„2 



=-1;- 4— log A 
ay 



we get 

v(^\y) = ^gb(2a + 3b)e'''y + 16/6^6^"^ + — . 
3 9 

Thus v^^^ can be obtained from v by the substitution a —a and g —g. It mean that 
v{y) is an SP-potential. 

B. Odd superpotential. In this case we have 

log^(y) = -ay - by^. 

A calculation yields 

v{y) = 4:{2by + af - ib. 

It is nothing but the harmonic oscillator and, of course an SP-potential: v^^^ — v + const. 

C. Exponential potential. Superpotential is W{(p) — 26e~'^, where 0(y) = log(a — by). 
Therefore 

It is a well known SP potential v^^^ — const x v. 



5 On models with cosmological expansion 

If 7^ in (2) then we have more complex system instead of (4) (in the case of the single 
brane at y = 0): 

Introducing t/j — A^, one can still reduce this system to the Schrodinger equation 

i^^{U{y)-pS{y)-X)^l^, 

where 

The reader may think that it is a very naive step, however this is not the case! In fact, for 
any given U{y) one can solve (5) to find ip and find V from (22). It is clear that solving 
(5) we should obtain a positive ip in order to find V, but it is not a problem, as has been 
demonstrated in sections 3 and 4. Indeed, one can construct a positive solution via the 
DT starting out from any simple solvable model, such as the RS model for example. 

But here wc face with another problem. To elucidate it, let us start out with an "RS 
potential" U = fi^, and consider the case A = 0. It is clear that the solution has the form 
(14). Using it as the prop function in (7), one gets the same metric A^.^ from section 3 
and, of course, with another potential. The main problem now is the kinetic energy of 
the scalar field, having assumed the form 

(0i'V 3^2^+6+ 3H^ I 

^ 2\a^ei'y ^b^e-i'y^ ^ 

It is clear that for large enough y, one gets 0^ < 0. On the other hand, if (see (16) 

(E^ sin^ 2q;± 



then on the branc and near it, there is no trouble. Therefore a simple way to avoid this 
problem is to deal with the finite volume, which is less then the value of y resulting in 
the negative kinetic term. Then one can use the DT and produce exact solutions in the 
brane world with cosmological expansion as well. 

6 Conclusion 

This approach can be easy generalized to the case of X-branes. To do it, we need to 
replace the term (35{y) in (5) by 

b 

and replace the two jump conditions (9) by 2K jump conditions at y = y^- This is 
technical matter, and the detail has been omitted. 

Another kettle of fish is that in the case of general position one cannot reconstruct the 
form of the potential V^(0), as a function of 0. Nevertheless, all these solutions will have 
the RS asymptotics, if constructed as has been described in section 3. The reader can 
easily obtain a great many exact solutions with the RS asymptotics, using the formulae 
from this section. ^ The DT gives a link between solvable problems, and one finds that 
most of (if not all) the exactly solvable potentials, from the harmonic oscillator to the 
finite-gap potentials [10], can be obtained via these transformations. The physical sense 
of these potentials in the brane world is not clear. Our main purpose here has been merely 
to demonstrate and advertise the DT as a powerful tool to manufacture exactly soluble 
potentials in the 5-D gravity with the bulk scalar field. It has been demonstrated that 
this method can be useful both in models without and with the cosmological expansion, 
although in the latter case the situation is more complex and not quite clear. 

There is notably one imperfection of this method: the DT would work for the 5-D 
gravity only. If > 5 then the warp factor is a multi-variable function. Unfortunately, 
there is no cogent and effective general DT theory in several dimensions [11]. 

But if D = 5 (and maybe if = 0) the DT is probably the best way to study the 
exactly solvable potentials. 
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